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Abstract. - From relaxation simulations of decoration-based quasicrystal structure models us- 
ing microscopically based interatomic pair potentials, we have calculated the (usually neglected) 
phonon-phason coupling constant. Its sign is opposite for the two alloys studied, i-AIMn and 
i-(Al,Cu)Li; a dimensionless measure of its magnitude relative to the phonon and phason elastic 
constants is of order 1/10, suggesting its effects are small but detectable. We also give a criterion 
for when phonon-phason effects are noticeable in diffuse tails of Bragg peaks. 



The elastic description of an icosahedral quasicrystal includes both "phonon" (ordinary) 
strain and "phason" strain. "Phason" strain parametrizes the way in which the local tile 
distribution deviates from icosahedral symmetry, as forced by the constraint of packing with 
surrounding tiles. Correspondingly there are two easily measured "phonon" elastic constants 
(ordinary bulk and shear modulus), and two "phason" elastic constants. The latter have been 
computed theoretically for various tilings pi 0, [|, || and recently measured experimentally || . 

There remains a cross-term with one elastic constant, the "phonon-phason" coupling r"^, 
which has not previously been measured experimentally or theoretically. The constant r" ,J - 
enters the shape of the diffuse tails around Bragg peaks |J, can drive a quasicrystal phase 
unstable with respect to a related crystal phase [[7| || , and affects the strain fields of dislocations 
in quasicrystals, which have been studied by diffraction contrast in transmission electron 
microscopy ||, |l0|]. Previous authors || |[ |l^, simply assumed T^' 1 - was small compared 
with the other elastic constants. 

The present work is the first attempt to compute r"'^ theoretically. As explained below, 
we extract it numerically from simulations in which we relax the atomic positions of a finite 
"periodic approximant" of the quasicrystal and measure its spontaneous shear distortion in 
response to the its intrinsic phason strain, for several different approximants. This was 
carried out for moderately realistic models of both z-AIMn and «-(Al,Cu)Li (modeled as a 
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pseudo-binary), representing the two major classes of quasicrystal. (Here "model" denotes a 
combination of a tiling, an atomic decoration rule, and a microscopically based pair potential.) 
In our conclusion, we crudely estimate ratios which determine the experimental observability 

off"'- 1 . 



Elastic theory. - This is formulated in terms of two kinds of strain fields: the "phonon" 



(ordinary) strain tensor has components 



(diUj + djUi)/2 where u(r) is the phonon 



displacement. The phason strain components are i?y = diVj. Here v(r) is the phason (also 
called "perp") displacement field, defined for quasicrystals and/or tilings; it can be constructed 
for any configuration of the tiles JlJ, and parameterizes the local imbalance in tiles 
of different orientations. Periodic approximant structures fl5f| of the quasicrystalline state 
necessarily acquire a non-zero phason strain. 

Group theory |^|, |l4| permits the following terms in the elastic free energy (in the notation 
of |Q, and adopting the summation convention): 



F 

Fphonon 
Fphason 
■phonon— phason 



Fphonon H~~ F'phason H~ ^phonon— phason 



(i) 



d s r- 



KiVijVij + K 2 \ v\ k - -VijUij + [(tvi2 + t 1 v 2 i) 2 + perms] 



d 3 rrii- L {[« u («n 



tu 2 2 + t 1 u 33 ) + 2u 23 (r 1 v 2 3 - tv 32 )] + perms} 



Here "perms" represents two terms obtained from the preceding term (within the same paren- 
thesis) by cyclic permutations of the indices (123). The term F p honon is the free energy of 
any isotropic harmonic solid with strain. It has two independent Lame constants A and /i, the 
latter being the shear modulus. Its "phason" analog is the F p h a son term with elastic constants 
K\ and K 2- The F p honon-phason term couples the shear strain with one representation of the 
phason strain. 

As in ordinary crystals, the elastic free energy describes the effects of thermal fluctuations 
on diffraction peaks, as well as possible soft-mode instabilities of the quasicrystalline state. 

One application of elastic theory is to dislocations in a quasicrystal, since any Burgers 
vector necessarily has both phonon and phason components. The smaller the ratio K\/fj,, the 
smaller the (expected) phonon component and the larger the phason component for a stable 
dislocation || [To). This might be modified by a nonzero T^ :± . 

Method for measuring phonon-phason coupling r"'^. 

Given any structure model, r"^ may be measured from the intrinsic distortion common 
to all approximant tilings. Its origin is the (known) nonzero uniform phason strain (vij) 
characteristic of each periodic crystalline approximant. Given (vij), the elastic theory then 
predicts a corresponding nonzero uniform (phonon) shear strain. Concretely, this means each 
tile spontaneously distorts a bit proportional to r" ,J - times (vij), so that tile edges in different 
directions are no longer quite related by icosahedral symmetry. 

Consider in particular an orthorhombic (pseudo-tetragonal) approximant with v\\ — V22 7^ 
U33, and all other Ujj = 0. We allow it to distort with fixed-volume shears Un = M22 = 
— e/2,U33 = e, so that the bulk modulus A does not contribute to the energy. Under this 
constraint, the minimum free energy occurs at 



2)i 



(«33 - Ull) 



(2) 
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Given a particular approximant, we can numerically compute the relaxed structural energy 
difference E re i(e) for several choices of shear magnitude e. Obtaining e m i n and /i by fitting 

E iel(e) = - fmm) 2 + const, (3) 
we finally deduce T"' 1 - using (||). 



Specification of the models. Even a crude estimate of the phonon-phason coupling 
requires sensible atomic models and realistic interatomic potentials jl6], [l7|], otherwise the 
structure is not very stable and likely to rearrange irreversibly after small perturbations. 
There are two important and different icosahedral structure types: the Frank-Kasper class 
and the Al-transition metal class, represented by i-(Al,Cu)Li and i-AlMn. We begin with the 
specification of the structural model, the first kind of input to our calculation. 

The atomic positions in the unit cell to be relaxed are determined by three partly indepen- 
dent structural factors: tiling, decoration, and approximant. Jig ]. The simplest icosahedral 
tiling geometry is the well-known three-dimensional rhombohedral tiling (3DRT) with edges 
along 5-fold symmetry axes |H| [l8|]. The Canonical Cell Tiling (CCT) [|l9| is an alternative 
tiling of 4 larger cells, which permits adjustments of the decoration to eliminate high-energy 
environments. 

In the AIMn system we have used one a representative decoration rule from each of three 



families of decoration rules for the CCT |16, p_7| These are labelled DF, DS, and LS in the 



scheme of |L6|, according to their tendency to dense (D) or loose (L) space- filling by the atoms, 
and their space-group symmetry (S=Simple Icosahedral or F=Face-Centered Icosahedral). Fj 
In the (Al,Cu)Li system, we only use the "Henley-Elser" (HE) decoration of the 3DRT |18) 
the version in which rhombohedra are grouped into disjoint rhombic dodecahedra wherever 
possible) . The HE decoration may also be applied to the CCT (which can be decomposed into 
a special case of the 3DRT [|i"9|). 

Having chosen the tiling type, we may select among the available anisotropic approximants. 
The lattice parameters for the unit cell of quasicrystal approximants are constrained to have 
only lengths which are powers of r times the quasi-lattice parameter [ fl5| . We simplify the 
conventional notation u p/q" Jig] for a rational approximant of r by only using u p" in our 
labelling. Thus, the anisotropic approximants are, ordered by increasing system size, 112, 221, 
223, 331, 332, ... §, f|J. In the CCT case, some of these cannot be tiled by the 4 CCT cells; 
in particular, there is no approximant between 331 and 335 ||. 

Phason strain will be measured in units where the the rhombohedron edge length an is taken 
to be unity in real space, as well as the corresponding edge in "perpendicular" space [p"3| . Thus 
v aa = t~ 3 , — t~ 4 , t~ 5 when the a entry in the approximant's label is 1,2, and 3 respectively. 

The second input in the specification of a model is the choice of pair potential. This is the 
only feasible approach to the energetics, since the system size of atomic structures in typical 
quasicrystal approximants is orders of magnitude larger than the maximum allowed for ab 
initio calculations. 

We adopted available well-tested potentials, constructed using linear response theory [ pl[ , 
for the binary AIMn Il7]] and the pseudo-binary (Al,Cu)Li (22| systems. The AIMn pair 



( 1 ) We used variants DF1.1, DS1.1, and LS1.1; we have not tried the other variants possib e 
within each family, which differ by changes of the chemistry or density on particular sites M, p7|; 
a representative LF1.1 from the fourth family of decorations in jl6| did not even show a quadratic 
E re \(t) function. 

( 2 ) Neither of these alloys is an equilibrium quasicrystal; the potentials for better quasicrystals 
were not yet available. 
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Table I. - Computed optimal shear e min , with shear modulus u and phonon-phason coupling r"' x in 
eV/atom. Text explains first three columns. 



Potential 


Structure 


Approximant 


^min 


M 


rii.-L 


AIMn 


CCT(DF) 


3 3 1 


0.03% 


6.7 


0.021 




CCT(DS) 


3 3 1 


0.11% 


6.3 


0.066 




CCT(LS) 


3 3 1 


0.22% 


7.6 


0.168 


(Al.Cu)Li 


CCT(HE) 


3 3 1 


-1.00% 


1.9 


-0.187 




3DRT(HE) 


1 1 2 


1.13% 


2.1 


-0.146 






2 2 1 


-0.81% 


2.2 


-0.110 






2 2 3 


-0.16% 


2.4 


-0.063 






3 3 1 


-0.80% 


2.2 


-0.177 






3 3 2 


0.16% 


1.8 


-0.047 



potentials involve the conventional perturbation theory for the nearly-free electrons in the 
Al-Al interaction, and Green's function treatment of Mn d-bands p3[ . For the (Al,Cu)Li 
system, ternary potentials were first constructed from pseudopotential-based perturbation 
theory and a tight-binding approximation of the Cu d-orbitals, and then the pseudo-binary 
potentials were obtained through proper averaging to reflect the chemical disorder in the 
system 0]. We call such pair potentials "realistically oscillating" JlTj because they feature 
Friedel oscillations out to several Angstroms. These potentials were tested in previous studies 
(see |l7j and J2^|); as a further check we applied the (Al,Cu)Li potentials to various crystal 
structures related to the Frank-Kasper family, and reproduced a sensible phase diagram 

Calculations and results. - We first constructed the atomic structure for each combination 
of model and pseudo-tetragonal approximant shown in table I. (Its first three columns identify 
the alloy and potential, the tiling type (and decoration), and the approximant, using notations 
explained above. The next larger approximant would have be "335" which was too large for 
our relaxations. 

For our selected pair potential, we always chose a cutoff slightly beyond 6A, so as to 
sufficiently include the first two potential wells |L7|. We applied to the cell a uniform fixed- 
volume shear e and then relaxed the positions (keeping the new cell shape fixed.) A modified 
conjugate gradient energy relaxation for the system converges to within 10~ 6 eV/atom. We 
sampled a grid of e values to derive the relaxed energy curve E rc i(e) which was fitted to (ph. 
Combining this with (|^) we extracted the shear modulus [i and phonon-phason coupling 
expressed in table I in units of eV/atom. 

For a given interatomic potential, the ordinary shear stiffness has reasonably consistent 
values. From table I, we first note /z « 7 eV/atom in i-AIMn or 2 eV/atom for i-(Al,Cu)Li. ( 3 ) 
The number density is about 6.5 atoms/aij 3 in i-AlMnJl(| or 7.8 atoms/a^ 3 in i-(Al,Cu)Li[|l8|. 
with aji about 4.6A or 5.05A, respectively. This gives the predictions /i w 75 GPa in i-AIMn 
or 20 GPa in i-(Al,Cu)Li, to be checked against the experimental values /z = 65 GPa in 
i-AlPdMn H or 41 GPa in i-AlCuLi §f|. 

Less surely, the table shows Ir'l-^l ~ 0.1 eV/atom. Although the magnitude of the phonon- 
phason coupling is ill-determined in table I, its sign is quite clear and not universal: it is positive 
for AIMn and negative for (Al,Cu)Li. We have no microscopic explanation for these signs, and 
so cannot tell whether they hold generally for Al-TM and Frank-Kasper type quasicrystals. 



( 3 ) Note that i-AIMn has a larger /i because its Al-Al potential has an energy scale about twice 
as big as in the (Al,Cu)Li case; in turn, this is due to the different electron densities. 
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The fitted rll ,J - for i-(Al,Cu)Li from different approximants varies by a factor ~ 4. We 
attribute this to finite-size effects due to the smallness of the tractable approximants: the 
strain relaxation might be especially sensitive to certain rare environments, occurring only one 
or two (or no) times in various approximants, producing big fluctuations in the inferred r" ,_L . 

We also relaxed all the 3DRT approximants of the (Al,Cu)Li structure using Roth's single- 
well Lennard- Jones (LJ) potential with non-additive radii |2g| ], rescaled to roughly match the 
first well of the realistic potentials. When the simulated binary alloy freezes under these LJ 
potentials, it is known to form a disordered quasicrystal of the Frank-Kasper class j28|. We 
found that T"' 1 - is much smaller in this case; our LJ results are tabulated in p7[ . 

Our calculation depended on the assumption that the optimal shear is a function only of 
the approximant's phason strain, and insensitive to tile rearrangements (which do not change 
phason strain) . Unfortunately, the only such comparison in our data set is between the 3DRT 
and CCT tilings in the (Al,Cu)Li 331 approximant. In that case (see Table 1), similar elastic 
constants were indeed fitted. 

We also adopted a crude and implausible "frozen" tiling assumption, that the energetic 
and entropic contributions to the elastic free energy are independent and linearly additive. 
The value of r" ,J ~ might be substantially renormalized (and be temperature-dependent) if we 
allowed atom relaxations and tile reshufflings simultaneously. 

Discussion. - We may now assess whether the phonon-phason term is strong enough to be 
measured experimentally, and even perhaps to affect the quasicrystal's stability. Near a Bragg 
vector G", the diffuse intensity (following Ishii's notation|8|) is proportional to 

(G ± ,B(q)- 1 G i )-2(Gll,A(q)- 1 C(q)B(q)- 1 G ± ) + (Gll,A(q)- 1 Gll) (4) 

We wish to estimate the relative importance of these three terms, representing respectively 
contributions from ordinary "phonons", from "phason" fluctuations, and from the phonon- 
phason cross term. The 3x3 matrices A(q), B(q), and C(q) are derived from the elastic 
free energy ([j]) and their formulas are given in eqs. (2.3)-(2.5) of ||]. Each matrix is a sum of 
terms proportional to the corresponding elastic constants and homogeneous of second order 
in |q|. Crudely speaking, the matrices scale as ||A(q)|| ~ /Li|q| 2 , ||B(q)|| ~ i^lqj 2 , and 
||C(q)|| ~ rll'^lqj 2 . (This scaling neglects their strong angular dependences on q/|q|; it also 
ignores the distinctions between the two phonon or the two phason elastic constants, since 
indeed A ~ 2/i in a typical solid, and \K 2 \/Ki = 0(1) for all cases quoted below.) 

We may conveniently express important criteria using ratios which all have the same 
dimensions as phason strain. First, the phason fluctuations - the first term in (|J) - surely 
dominate that sum, as assumed in the analysis of [^J, thus 

|G^|/|Gll|>|rH^|/ M . (5) 

Second, for the phonon-phason coupling to be measurable, it ought to be the dominant 
correction, i.e. the second term in (|j) is the next largest. Thus 

|G^/|Gll|>AV|r" a |- (6) 



Adopting the simplest sort of random-tiling scenario p2|, all of the phason elastic term 

-3 

it 



is part of the tiling's configurational entropy. Simulations find K\/T = 0.81 a R 3 with 



K 2 /K 2 = 0.61 for the 3DRT g, §. For a model approximating the CCT [|, they find K 1 /T 
is (at most) 2.5a^ 3 with K2/K1 ~ —0.75. ( 4 ) Furthermore, phason fluctuations fall out of 



( 4 ) A less reliable calculation, but for the true CCT, found K\ w a R 
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equilibrium not so far below the melting temperature T m , so we will replace T — * T rn m 10 3 K 
in these estimates. Thus, using the figures mentioned earlier for an and the number density, a 
crude estimate for either alloy is K\ ~ 0.03 eV/atom, in the units of table I. (The experimental 
value of K\ in i-AlPdMn j2{| seems to be a factor of 10 smaller.) 

The above quoted estimates for fj, (for i-AlMn), r"^, and K\, when inserted into the 
right-hand sides of (||) and (^), give respectively 0.015 and 0.3. Thus a Bragg peak may be 



found such that the T^' 1 - term is the first correction. This correction would cause differences 
in shape of the (otherwise similar) diffuse scattering contours of two Bragg peaks with G" in 
the same direction, as observed recently |^9| . 

There is an additional correction (see e.g. (gj) to every term of (Q), scaling with the relative 
order 

r = irll'^/Xi/i (7) 

(provided we assume K% — 0{K\) and A = O(p), as in the previous discussion). Evidently r 
is the (dimensionless) ratio of the right hand sides of (|J) and (||). Using the elastic constants 
suggested above for i-AIMn and i-(Al,Cu)Li gives r ~ 0.05 and r ~ 0.2, respectively, uncertain 
by an order of magnitude either way. This correction provides a second way that r" ,J ~ affects 
the diffuse scattering. Thus too, a version of r shows up in the criterion for nondivergence of 
the fluctuations (and diffuse scattering), Widom's eq. (7): C T (rH a ) 2 /{\+2^)(K^K 2 /Z) < 1 
where C_ = 4/9 and C+ = 4. 

The same ratio r also determines the importance of T"'' L in the onset of instability to a 
uniform phason strain, as discussed in jjj and [B. Such phason-elastic instabilities might 
generate the approximant phases - e.g. rhombohedral R(AlCuFe) - typically occurring for 
compositions close to the quasicrystal's |@]. The stability criterion is given in Widom's eq. 
(5) @: 3|rll' ± | 2 /(Ki - AK 2 /3)n < 1- If K2/K1 « 0.6, as in the random 3DRT § or in 
i-AlCuFe [^9|, the criterion says r < 0.07. If our estimate of r is valid for real i-AlCuFe, then 
it is close to a phonon-phason driven instability. 

The smallness of the phonon-phason constant has been rationalized by the notion that it 
measures the "frustration" of a material, defined as the mismatch of the ideal interatomic 
spacings and the well radii of the potentials |0|. In the unfrustrated case, different terms 
of the interactions are independently satisfied by the atoms, which occupy tiling-like special 
positions. This permits good quasicrystal order, since the tile-tile energy may be low for more 
than one way of packing tiles along a face. At the same time, the lack of unbalanced forces on 
the atoms means there is less tendency for tile shapes to distort depending on the tile packing. 
In this way the "binary tiling" , a two-dimensional atomic toy model [j30| , has a quite small 

However, our results appear to contradict the above argument: among the three AIMn deco- 
rations of the CCT 331 approximant, the "LS" decoration is best adapted to the structure [ p"7j |, 
yet it shows the highest Til'- 1 . (We suppose this is somehow due to the small voids common 
in the "LS" structure, which allow greater atomic rearrangement in response to phonon or 
phason shears.) Similarly, on the 3DRT tilings with the i-(Al,Cu)Li structure, the i-(Al,Cu)Li 
potentials in Table I give a much larger phonon-phason coupling than the poorer LJ model 
(computed in p7|). 

In conclusion, we have numerically estimated the phonon-phason coupling r" ,J - for repre- 
sentatives of both major classes of icosahedral quasicrystal. We found it had a non-universal 
sign, and (with very large uncertainties) it was large enough in magnitude to be observable 

( B ) Bragg peaks exist with arbitrarily small |G"|/|G X |, but with exponentially small intensi- 
ties g. 
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and perhaps even to dominate the mechanism of the quasicrystal's instability to a uniform 
phason strain. To improve the calculations, a stable quasicrystal should be modeled and the 
fluctuations in small samples (due to variations in tiling configurations and in decoration rules) 
must be better understood. 

*** 
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